Role of SARS-CoV-2 mutations in the evolution of the COVID-19 pandemic.

OBJECTIVES: The SARS-CoV-2 pandemic and large-scale genomic surveillance provided an exceptional opportunity to analyze mutations that appeared over three years in viral genomes. Here we studied mutations and their epidemic consequences for SARS-CoV-2 genomes from our center. METHODS: We analyzed 61,397 SARS-CoV-2 genomes we sequenced from respiratory samples for genomic surveillance. Mutations frequencies were calculated using Nextclade, Microsoft Excel, and an in-house Python script. RESULTS: A total of 22,225 nucleotide mutations were identified, 220 (1.0%) being each at the root of ≥836 genomes, classifying mutations as 'hyperfertile'. Two seeded the European pandemic: P323L in RNA polymerase, associated with an increased mutation rate, and D614G in spike that improved fitness. Most 'hyperfertile' mutations occurred in areas not predicted with increased virulence. Their mean number was 8±6 (0-22) per 1000 nucleotides per gene. They were 3.7-times more frequent in accessory than informational genes (13.8 versus 3.7/1000 nucleotides). Particularly, they were 4.1-times more frequent in ORF8 than in the RNA polymerase gene. Interestingly, stop codons were present in 97 positions, almost only in accessory genes, including ORF8 (21/100 codons). CONCLUSIONS: most 'hyperfertile' mutations did not predict emergence of a new epidemic, and some were stop codons indicating the existence of so-named 'non-virulence' genes.

Resistance-associated mutations to the anti-SARS-CoV-2 agent nirmatrelvir: Selection not induction.

Mutations associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resistance to antiprotease nirmatrelvir were reported. We aimed to detect them in SARS-CoV-2 genomes and quasispecies retrieved in our institute before drug availability in January 2022 and to analyze the impact of mutations on protease (3CLpro) structure. We sought for 38 3CLpro nirmatrelvir resistance mutations in a set of 62 673 SARS-CoV-2 genomes obtained in our institute from respiratory samples collected between 2020 and 2023 and for these mutations in SARS-CoV-2 quasispecies for 90 samples collected in 2020, using Python. SARS-CoV-2 protease with major mutation E166V was generated with Swiss Pdb Viewer and Molegro Molecular Viewer. We detected 22 (58%) of the resistance-associated mutations in 417 (0.67%) of the genomes analyzed; 325 (78%) of these genomes had been obtained from samples collected in 2020-2021. APOBEC signatures were found for 12/22 mutations. We also detected among viral quasispecies from 90 samples some minority reads harboring any of 15 nirmatrelvir resistance mutations, including E166V. Also, we predicted that E166V has a very limited effect on 3CLpro structure but may prevent drug attachment. Thus, we evidenced that mutations associated with nirmatrelvir resistance pre-existed in SARS-CoV-2 before drug availability. These findings further warrant SARS-CoV-2 genomic surveillance and SARS-CoV-2 quasispecies characterization.

The return of the "Mistigri" (virus adaptative gain by gene loss) through the SARS-CoV-2 XBB.1.5 chimera that predominated in 2023.

Severe acute respiratory syndrome coronavirus 2 XBB.1.5 is the first recombinant lineage to predominate at the country and global scales. Very interestingly, like the Marseille-4B subvariant (or B.1.160) and the pandemic variant B.1.1.7 (or Alpha) previously, it has its ORF8 gene inactivated by a stop codon. We aimed here to study the distribution of stop codons in ORF8 of XBB.1.5 and non-XBB.1.5 genomes. We identified that a stop codon was present at 89 (74%) ORF8 codons in ≥1 of 15 222 404 genomes available in GISAID. The mean proportion of genomes with a stop codon per codon was 0.11% (range, 0%-7.8%). In addition, a stop codon was detected at 15 (12%) codons in at least 1000 genomes. These 15 codons are notably located on seven stem-loop hairpin regions and in the signal peptide region for the case of the XBB.1.5 lineage (codon 8). Thus, it is very likely that stop codons in ORF8 gene contributed on at least three occasions and independently during the pandemic to the evolutionary success of a lineage that became transiently predominant. Such association of gene loss with evolutionary success, which suits the recently described Mistigri rule, is an important biological phenomenon very unknown in virology while largely described in cellular organisms.

Emergence of a second SARS-CoV-2 variant with a tremendous genetic leap from its ancestors.

The on-going emergence of the Omicron BA.2.86 variant is one of the major events in SARS-CoV-2 genetic evolution that remain enigmatic regarding the overall virus mutation rate, together with the emergence of the initial Omicron variant, BA.1. Indeed, the genomes of the Omicron BA.2.86 lineage, an offspring of the second major Omicron subvariant, BA.2, harbor 39 additional mutations in the spike compared to this ancestor. Here we comment on the phylogeny of BA.2.86, on the positions, and frequencies in other SARS-CoV-2, of mutations in its spike, and on the structural model of this protein that concentrates most of BA.2.86 additional mutations.

Sequencing of monkeypox virus from infected patients reveals viral genomes with APOBEC3-like editing, gene inactivation, and bacterial agents of skin superinfection.

A large outbreak of Monkeypox virus (MPXV) infections has arisen in May 2022 in nonendemic countries. Here, we performed DNA metagenomics using next-generation sequencing with Illumina or Nanopore technologies for clinical samples from MPXV-infected patients diagnosed between June and July 2022. Classification of the MPXV genomes and determination of their mutational patterns were performed using Nextclade. Twenty-five samples from 25 patients were studied. A MPXV genome was obtained for 18 patients, essentially from skin lesions and rectal swabbing. All 18 genomes were classified in clade IIb, lineage B.1, and we identified four B.1 sublineages (B.1.1, B.1.10, B.1.12, B.1.14). We detected a high number of mutations (range, 64-73) relatively to a 2018 Nigerian genome (genome GenBank Accession no. NC_063383.1), which were harbored by a large part of a set of 3184 MPXV genomes of lineage B.1 recovered from GenBank and Nextstrain; and we detected 35 mutations relatively to genome ON563414.3 (a B.1 lineage reference genome). Nonsynonymous mutations occurred in genes encoding central proteins, among which transcription factors and core and envelope proteins, and included two mutations that would truncate a RNA polymerase subunit and a phospholipase d-like protein, suggesting an alternative start codon and gene inactivation, respectively. A large majority (94%) of nucleotide substitutions were G > A or C > U, suggesting the action of human APOBEC3 enzymes. Finally, >1000 reads were identified as from Staphylococcus aureus and Streptococcus pyogenes for 3 and 6 samples, respectively. These findings warrant a close genomic monitoring of MPXV to get a better picture of the genetic micro-evolution and mutational patterns of this virus, and a close clinical monitoring of skin bacterial superinfection in monkeypox patients.

New Genotype of Coxiella burnetii Causing Epizootic Q Fever Outbreak in Rodents, Northern Senegal.

In Senegal, Coxiella burnetii, which causes Q fever, has often been identified in ticks and humans near livestock, which are considered to be reservoirs and main sources of infection. We describe the emergence of C. burnetii in rodents, not previously known to carry this pathogen, and describe 2 new genotypes.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule.

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

Pangenomic analysis of Coxiella burnetii unveils new traits in genome architecture.

Coxiella burnetii is the etiological agent of Q fever, a worldwide zoonosis able to cause large outbreaks. The disease is polymorphic. Symptomatic primary infection is named acute Q fever and is associated with hepatitis, pneumonia, fever, and auto-immune complications while persistent focalized infections, mainly endocarditis, and vascular infections, occur in a minority of patients but are potentially lethal. In order to evaluate the genomic features, genetic diversity, evolution, as well as genetic determinants of antibiotic resistance, pathogenicity, and ability to cause outbreaks of Q fever, we performed a pangenomic analysis and genomic comparison of 75 C. burnetii strains including 63 newly sequenced genomes. Our analysis demonstrated that C. burnetii has an open pangenome, unique genes being found in many strains. In addition, pathogenicity islands were detected in all genomes. In consequence C. burnetii has a high genomic plasticity, higher than that of other intracellular bacteria. The core- and pan-genomes are made of 1,211 and 4,501 genes, respectively (ratio 0.27). The core gene-based phylogenetic analysis matched that obtained from multi-spacer typing and the distribution of plasmid types. Genomic characteristics were associated to clinical and epidemiological features. Some genotypes were associated to specific clinical forms and countries. MST1 genotype strains were associated to acute Q fever. A significant association was also found between clinical forms and plasmids. Strains harboring the QpRS plasmid were never found in acute Q fever and were only associated to persistent focalized infections. The QpDV and QpH1 plasmids were associated to acute Q fever. In addition, the Guyanese strain CB175, the most virulent strain to date, exhibited a unique MST genotype, a distinct COG profile and an important variation in gene number that may explain its unique pathogenesis. Therefore, strain-specific factors play an important role in determining the epidemiological and clinical manifestations of Q fever alongside with host-specific factors (valvular and vascular defects notably).

Quasispecies Analysis of SARS-CoV-2 of 15 Different Lineages during the First Year of the Pandemic Prompts Scratching under the Surface of Consensus Genome Sequences.

The tremendous majority of SARS-CoV-2 genomic data so far neglected intra-host genetic diversity. Here, we studied SARS-CoV-2 quasispecies based on data generated by next-generation sequencing (NGS) of complete genomes. SARS-CoV-2 raw NGS data had been generated for nasopharyngeal samples collected between March 2020 and February 2021 by the Illumina technology on a MiSeq instrument, without prior PCR amplification. To analyze viral quasispecies, we designed and implemented an in-house Excel file ("QuasiS") that can characterize intra-sample nucleotide diversity along the genomes using data of the mapping of NGS reads. We compared intra-sample genetic diversity and global genetic diversity available from Nextstrain. Hierarchical clustering of all samples based on the intra-sample genetic diversity was performed and visualized with the Morpheus web application. NGS mapping data from 110 SARS-CoV-2-positive respiratory samples characterized by a mean depth of 169 NGS reads/nucleotide position and for which consensus genomes that had been obtained were classified into 15 viral lineages were analyzed. Mean intra-sample nucleotide diversity was 0.21 ± 0.65%, and 5357 positions (17.9%) exhibited significant (>4%) diversity, in ≥2 genomes for 1730 (5.8%) of them. ORF10, spike, and N genes had the highest number of positions exhibiting diversity (0.56%, 0.34%, and 0.24%, respectively). Nine hot spots of intra-sample diversity were identified in the SARS-CoV-2 NSP6, NSP12, ORF8, and N genes. Hierarchical clustering delineated a set of six genomes of different lineages characterized by 920 positions exhibiting intra-sample diversity. In addition, 118 nucleotide positions (0.4%) exhibited diversity at both intra- and inter-patient levels. Overall, the present study illustrates that the SARS-CoV-2 consensus genome sequences are only an incomplete and imperfect representation of the entire viral population infecting a patient, and that quasispecies analysis may allow deciphering more accurately the viral evolutionary pathways.

SARS-CoV-2 quasi-species analysis from patients with persistent nasopharyngeal shedding.

At the time of a new and unprecedented viral pandemic, many questions are being asked about the genomic evolution of SARS-CoV-2 and the emergence of different variants, leading to therapeutic and immune evasion and survival of this genetically highly labile RNA virus. The nasopharyngeal persistence of infectious virus beyond 17 days proves its constant interaction with the human immune system and increases the intra-individual mutational possibilities. We performed a prospective high-throughput sequencing study (ARTIC Nanopore) of SARS-CoV-2 from so-called "persistent" patients, comparing them with a non-persistent population, and analyzing the quasi-species present in a single sample at time t. Global intra-individual variability in persistent patients was found to be higher than in controls (mean 5.3%, Standard deviation 0.9 versus 4.6% SD 0.3, respectively, p < 0.001). In the detailed analysis, we found a greater difference between persistent and non-persistent patients with non-severe COVID 19, and between the two groups infected with clade 20A. Furthermore, we found minority N501Y and P681H mutation clouds in all patients, with no significant differences found both groups. The question of the SARS-CoV-2 viral variants' genesis remains to be further investigated, with the need to prevent new viral propagations and their consequences, and quasi-species analysis could be an important key to watch out.

First Detection of the SARS-CoV-2 Omicron BA.5/22B in Monaco.

The Omicron BA.5/22B variant has been designated as a "variant of concern" by the World Health Organization. We describe, here, the first evidence in Monaco of infection with an Omicron BA.5/22B variant, probably imported from the Republic of Seychelles, harboring a rare combination of non-BA.5/22B signature amino acid changes. SARS-CoV-2 neutralizing antibodies were measured with a surrogate virus neutralization test. SARS-CoV-2 genotype screening was performed on nasopharyngeal samples with a multiplex qPCR assay. The SARS-CoV-2 genome was obtained by next-generation sequencing with the Illumina COVID-seq protocol, then assembly using bioinformatics pipelines and software was performed. The BA.5/22B spike protein structure was obtained by molecular modeling. Two spouses were SARS-CoV-2-diagnosed the day they returned from a one-week trip in the Republic of Seychelles. SARS-CoV-2 qPCR screening for variant-specific mutations identified an Omicron variant BA.1/21K, BA.4/22A, or BA.5/22B. A SARS-Co-2 BA.5/22B variant genome was recovered from one of the spouses. Aside from BA.5/22B-defining amino acid substitutions, four other amino acid changes were encoded including Q556K in ORF1a, K2557R in ORF1b, and A67V and A829T in spike; only 13 genomes in sequence databases harbored these four mutations concurrently. Structural analysis of this BA.5/22B variant predicted that A829T in spike may result in a compaction that may affect conformational plasticity. Overall, our findings warrant performing genome-based genotypic surveillance to survey accurately the emergence and circulation of SARS-CoV-2 variants worldwide and point out that their first occurrence in a country is often through international travel despite implemented countermeasures.

A 21L/BA.2-21K/BA.1 "MixOmicron" SARS-CoV-2 hybrid undetected by qPCR that screen for variant in routine diagnosis.

Among the multiple SARS-CoV-2 variants identified since summer 2020, several have co-circulated, creating opportunities for coinfections and potentially genetic recombinations that are common in coronaviruses. Viral recombinants are indeed beginning to be reported more frequently. Here, we describe a new SARS-CoV-2 recombinant genome that is mostly that of a Omicron 21L/BA.2 variant but with a 3' tip originating from a Omicron 21K/BA.1 variant. Two such genomes were obtained in our institute from adults sampled in February 2022 in university hospitals of Marseille, southern France, by next-generation sequencing carried out with the Illumina or Nanopore technologies. The recombination site was located between nucleotides 26,858-27,382. In the two genomic assemblies, mean sequencing depth at mutation-harboring positions was 271 and 1362 reads and mean prevalence of the majoritary nucleotide was 99.3 ± 2.2% and 98.8 ± 1.6%, respectively. Phylogeny generated trees with slightly different topologies according to whether genomes analyzed were depleted or not of the 3' tip. This 3' terminal end brought in the Omicron 21L/BA.2 genome a short transposable element of 41 nucleotides named S2m that is present in most SARS-CoV-2 except a few variants among which the Omicron 21L/BA.2 variant and may be involved in virulence. Importantly, this recombinant is not detected by currently used qPCR that screen for variants in routine diagnosis. The present observation emphasizes the need to survey closely the genetic pathways of SARS-CoV-2 variability by whole genome sequencing, and it could contribute to gain a better understanding of factors that lead to observed differences between epidemic potentials of the different variants.

Culture and identification of a "Deltamicron" SARS-CoV-2 in a three cases cluster in southern France.

Multiple SARS-CoV-2 variants have successively, or concomitantly spread worldwide since the summer of 2020. A few co-infections with different variants were reported and genetic recombinations, common among coronaviruses, were reported or suspected based on co-detection of signature mutations of different variants in a given genome. Here we report three infections in southern France with a Delta 21J_AY.4-Omicron 21K/BA.1 "Deltamicron" recombinant. The hybrid genome harbors signature mutations of the two lineages, supported by a mean sequencing depth of 1163-1421 reads and a mean nucleotide diversity of 0.1%-0.6%. It is composed of the near full-length spike gene (from codons 156-179) of an Omicron 21K/BA.1 variant in a Delta 21J/AY.4 lineage backbone. Importantly, we cultured an isolate of this recombinant and sequenced its genome. It was observed by scanning electron microscopy. As it is misidentified with current variant screening quantitative polymerase chain reaction (qPCR), we designed and implemented for routine diagnosis a specific duplex qPCR. Finally, structural analysis of the recombinant spike suggested its hybrid content could optimize viral binding to the host cell membrane. These findings prompt further studies of the virological, epidemiological, and clinical features of this recombinant.

SARS-CoV-2 reinfection and COVID-19 severity.

SARS-CoV-2 reinfection rate is low. The relative severity of the first and second episodes of infection remains poorly studied. In this study, we aimed at assessing the frequency of SARS-CoV-2 reinfections and comparing the severity of the first and second episodes of infection. We retrospectively included patients with SARS-CoV-2 positive RT-PCR at least 90 days after clinical recovery from a COVID-19 episode and with at least one negative RT-PCR after the first infection. Whole genome sequencing and variant-specific RT-PCR were performed and clinical symptoms and severity of infection were retrospectively documented from medical files. A total of 209 COVID-19 reinfected patients were identified, accounting for 0.4% of positive cases diagnosed from 19 March 2020 to 24 August 2021. Serology was performed in 64 patients, of whom 39 (60.1%) had antibodies against SARS-CoV-2 when sampled at the early stage of their second infection. Only seven patients (3.4%) were infected twice with the same variant. We observed no differences in clinical presentation, hospitalization rate, and transfer to ICU when comparing the two episodes of infections. Our results suggest that the severity of the second episode of COVID-19 is in the same range as that of the first infection, including patients with antibodies.

First cases of infection with the 21L/BA.2 Omicron variant in Marseille, France.

The SARS-CoV-2 21K/BA.1, 21L/BA.2, and BA.3 Omicron variants have recently emerged worldwide. To date, the 21L/BA.2 Omicron variant has remained very minority globally but became predominant in Denmark instead of the 21K/BA.1 variant. Here, we describe the first cases diagnosed with this variant in south-eastern France. We identified 13 cases using variant-specific qPCR and next-generation sequencing between 28/11/2021 and 31/01/2022, the first two cases being diagnosed in travelers returning from Tanzania. Overall, viral genomes displayed a mean (±standard deviation) number of 65.9 ± 2.5 (range, 61-69) nucleotide substitutions and 31.0 ± 8.3 (27-50) nucleotide deletions, resulting in 49.6 ± 2.2 (45-52) amino acid substitutions (including 28 in the spike protein) and 12.4 ± 1.1 (12-15) amino acid deletions. Phylogeny showed the distribution in three different clusters of these genomes, which were most closely related to genomes from England and South Africa, from Singapore and Nepal, or from France and Denmark. Structural predictions highlighted a significant enlargement and flattening of the surface of the 21L/BA.2 N-terminal domain of the spike protein compared to that of the 21K/BA.1 Omicron variant, which may facilitate initial viral interactions with lipid rafts. Close surveillance is needed at global, country, and center scales to monitor the incidence and clinical outcome of the 21L/BA.2 Omicron variant.

Introduction of the SARS-CoV-2 Beta variant from Comoros into the Marseille geographical area.

BACKGROUND: We describe the epidemiology of the first cases diagnosed in our institute of infections with the SARS-CoV-2 Beta variant and how this variant was imported to Marseille. METHODS: The Beta variant was identified based on analyses of sequences of viral genomes or of a spike gene fragment obtained by next-generation sequencing using Illumina technology, or by a real-time reverse-transcription-PCR (qPCR) specific of the Beta variant. RESULTS: The first patient diagnosed as infected with the SARS-CoV-2 Beta variant was sampled on January 15, 2021. Twenty-nine patients were diagnosed in January 2021 (two weeks). Fifteen (52%) patients were of Comorian nationality. Eight (28%) had travelled abroad, including six who had returned from Comoros. Phylogeny based on SARS-CoV-2 genomes from 11 of these patients and their best BLAST hits from the GISAID database showed that seven patients, including the four returning from Comoros, were clustered with 27 other genomes from GISAID that included the six first Beta variant genomes described in Comoros in January 2021. CONCLUSIONS: Our analyses highlight that, as for the case of other SARS-CoV-2 variants that have been diagnosed in Marseille, the Beta variant was imported to Marseille through travel from abroad. It had limited spread in our geographical area.

Emergence in southern France of a new SARS-CoV-2 variant harbouring both N501Y and E484K substitutions in the spike protein.

SARS-CoV-2 variants have become a major virological, epidemiological, and clinical concern, particularly with regard to the risk of escape from vaccine-induced immunity. Here, we describe the emergence of a new variant, with the index case returning from travel in Cameroon. For 13 SARS-CoV-2-positive patients living in the same geographical area of southeastern France, a qPCR test for screening variant-associated mutations showed an atypical combination. The genome sequences were obtained by next-generation sequencing with Oxford Nanopore Technologies on GridION instruments within about 8 h. Analysis revealed 46 nucleotide substitutions and 37 deletions, resulting in 30 amino acid substitutions and 12 deletions. Fourteen of the amino acid substitutions, including N501Y and E484K, and nine deletions are located in the spike protein. This genotype pattern led to the establishment of a new Pangolin lineage, named B.1.640.2, that is a phylogenetic sister group to the old B.1.640 lineage, which has now been renamed B.1.640.1. The lineages differ by 25 nucleotide substitutions and 33 deletions. The combination of mutations in these isolates and their phylogenetic position indicate, based on our previous definition, that they represent a new variant, which we have named "IHU". These data are a further example of the unpredictability of the emergence of SARS-CoV-2 variants, and of their possible introduction into a given geographical area from abroad.

Limited spread of a rare spike E484K-harboring SARS-CoV-2 in Marseille, France.

We detected SARS-CoV-2 of PANGO lineage R.1 with the spike substitution E484K in three patients. Eleven other sequences in France and 8,831 worldwide were available from GISAID, 92% originating from Japan. The three genome sequences from our institute were phylogenetically closest to another from Guinea-Conakry, where one of the patients had travelled. These viruses did not exhibit any unusual features in cell culture. Spike structural predictions indicated a 1.3-time higher transmissibility index than for the globally spread B.1.1.7 variant but also an affinity loss for gangliosides that might have slowed dissemination. The spread of new SARS-CoV-2 mutants/variants is still not well understood and therefore difficult to predict, and this hinders implementation of effective preventive measures, including adapted vaccines.

Occurrence of a substitution or deletion of SARS-CoV-2 spike amino acid 677 in various lineages in Marseille, France.

Great concerns have been raised about SARS-CoV-2 variants over the past six months. At the end of 2020, an increasing incidence of spike substitutions Q677H/P was described in the USA, which involved six independent lineages. We searched for changes to this amino acid in the sequence database of SARS-CoV-2 genomes obtained at the IHU Méditerranée Infection (Marseille, France) from 3634 patients sampled between February 2020 and April 2021. In seven genomes (0.2%), we found a deletion of five amino acids at spike positions 675-679 (QTQTN) including Q677, and in 76 genomes (2.3%) we found a Q677H substitution. The 83 genomes were classified in ten different Pangolin lineages. Genomes with a spike Q677 deletion were obtained from respiratory samples collected in six cases between 28 March 2020 and 12 October 2020 and in one case on 1 February 2021. The Q677H substitution was found in genomes all obtained from respiratory samples collected from 19 January 2021 and were classified in seven different lineages. Most of these genomes (41 cases) were of UK variant. Two others were classified in the B.1.160 Pangolin lineage (Marseille-4 variant) which was first detected in July 2020 in our institute but was devoid of this substitution until 19 January 2021. Also, eight genomes were classified in the A.27/Marseille-501 lineage which was first detected in our institute in January 2021 and which either harboured or did not harbour the Q677H substitution. Thus, the spike Q677H substitution should be considered as another example of convergent evolution, as it is the case of spike substitutions L18F, E484K, L452R, and N501Y which also independently appeared in various lineages.